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The differential branching ratio, forward-backward asymmetry, CP asym- 
metry and lepton polarization for a B-meson to decay to strange hadronic final 
states and a r+r" pair in a CP softly broken two Higgs doublet model are 
computed. It is shown that contributions of neutral Higgs bosons to the decay 
are quite significant when tan /? is large. And it is proposed to measure the 
direct CP asymmetry in back-forward asymmetry. 
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I. INTRODUCTION 



The origin of the CP violation has been one of main issues in high energy physics since 
the discovery of the CP violation in the Kq — Kq systerm in 1964 [|I|. The measurements of 
electric dipole moments of the neutron and electron and the matter-antimatter asymmetry 
in the universe indicate that one needs new sources of CP violation in addition to the 
CP violation come from CKM matrix, which has been one of motivations to search new 
theoretical models beyond the standard model (SM). 

The minimal extension of the SM is to enlarge the Higgs sectors of the SM [Q. It has been 
shown that if one adheres to the natural flavor conservation (NFC) in the Higgs sector, then 
a minimum of three Higgs doublets are necessary in order to have spontaneous CP violations 
0. However, the constraint can be evaded if one allows the real and image parts of 05*' 02 
have different self-couplings and adds a linear term of 'Re{(f)i (j)2) in the Higgs potential (see 
below Eq. (^ with m\ = 0). Then, one can construct a CP spontaneously broken two 
Higgs doublet (2HDM), which is the minimal and the most "economical" one |^ among the 
extensions of the SM that provide new source of CP violation. Furthermore, in addition to 
the above terms, if one adds a linear term of Im(0j''02), then one has a CP softly broken 
2HDM §. 

Flavor changing neutral current (FCNC) transitions B Xg'j and B —>■ Xgl^l^ provide 
testing grounds for the SM at the loop level and sensitivity to new physics. Rare decays 
B Xgl^l^il = e, fi) have been extensively investigated in both SM and the beyond [^0. 
In these processes contributions from exchanging neutral Higgs bosons (NHB) can be safely 
neglected because of smallness of = e,/i). The inclusive decay B X^r+r^ has 

also been investigated in the SM, the model II 2HDM and SUSY models with and without 

* Comparing the Model III 2HDM in which CP is explicitly violated, the CP spontaneously 
broken 2HDM has only two new parameters besides the masses of the Higgs bosons in the large 
tan (3 limit (see below) . In this sense it is the most " economical" . 



1 



including the contributions of NHB [p|-p!l[]. In this note we investigate the inclusive decay 
B XgT^r^ with emphasis on CP violation effect in a CP softly broken 2HDM, which 
we shall call Model IV hereafter for the sake of simplisity. We consider the Model IV in 
which the up-type quarks get masses from Yukawa couplings to the one Higgs doublet H2 
and down-type quarks and leptons get masses from Yukawa couplings to the another Higgs 
doublet Hi. The Higgs boson couplings to down- type quarks and leptons depend on only 
the CP violated phase ^ which comes from the expectation value of Higgs and the ratio 
tg(3 = ^ in the large tg(3 limit (see next section), which are the free parameters in the 
model. Because the couplings of the charged Higgs to fermions in Model IV are the same as 
those in the model II, the constraints on tan (3 due to effects arising from the charged Higgs 
are the same as those in the model II. Constraints on tg(3 from K — K and B — B mixing, 



r(6 — > S7),r(6 ctUt) and Ri, have been given ||T2 



0.7<<«/J<0.52(^) (1) 



(and the lower limit ttih^ > 2QQGev has also been given in the ref. W^)- It is obvious 



that the contributions from exchanging neutral Higgs bosons now is enhanced roughly by a 
factor of tg'^(3 and can compete with those from exchanging 7, Z when tg(3 is large enough. 
Because the CP violation effects in 5 — > X^r+r^ come from the couplings of NHB to leptons 
and quarks, we shall be interested in the large tan (3 limit in this note. The constraints on 
^ can be obtained from the electric dipole moments (EDM) of the neutron and electron, 
which will be analysed in the next section. 

II. MODEL DESCRIPTION 

Consider two complex y = I, SU {2)^ doublet scalar fields, 0i and 02- The Higgs potential 
which spontaneously breaks SU(2) x f/(l) down to U{1)em can be written in the following 
form [|[: 

V^(01, 02) = E H<Pt<Pr + H^t<t>if] 
i=l,2 



+ A3[(0^0l)(0^02)] + A4[(0+02)(0^0l)] 

+A5[Re(0+02)]' + A6[Im(0+02)]' 



(2) 



Hermiticity requires that all parameters are real. The potential is CP softly broken due to 
the presence of the term m\Im{(j)i(j)2)- It is easy to see that the minimum of the potential 
is at 



< 01 >= 



' ^ 



< 02 > = 



V2e 



(3) 



thus breaking SU{2) x f/(l) down to U{1)em and simutaneously breaking CP, as desired. 
It should be noticed that only for A5 7^ Ag, the phase ^ can't rotated away as usual, which 
breaks the CP-conservation. If m|=0 in (2) then the potential is CP invarint. It has been 
shown that the CP spontaneously breaking happens at (3) ||13|. We limit ourself to the case 
of m\ 7^ in the paper and shall investigate the m\=Q case in a separate paper |T^. 

In the following we will work out the mass spectrum of the Higgs boson. For charged 
components, the mass-squared matrix for negtive states is 



A. 



V 



-ViV2e 



(4) 



Diagonalizing the mass-squared matrix results in one zero-mass Goldstone state: 



G = e'^ sin /?02 + cos ^(pi 



(5) 



and one massive charged Higgs boson state: 



H 



e*^ cos / 



niH- 



sm , 
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\Xi\{vl + vi), 



(6) 
(7) 



where tan /5 = V2/V1. Correspondingly we could also get the positive states G'^ and with 
the same masses zero and l^iUvf + V2), respectively. 



For neutral Higgs components, because CP-conservation is breaking, the mass-squared 
matrix is 4 x 4, which could not be simply separated into two 2x2 matrices as usual. 
However, in the case of large tan j3 which is we intrested in, the neutral parts can be written 
as separately two 2x2 matrices and one of them is 



/a. 



v 



A5+A6 + (A6-A5)cos(2g) 
2 

(A6-A5)sin(2g) 
2 



(A6~A5)sm(2g) 
2 

A5+A6+(A5~A6)cos(2g) 
2 



(8) 



Diagonalizing the Higgs boson mass-squared matrix results in two eigenstates: 

-V2 



Re0? 



(9) 



with masses 



m^o = XqvI, (10) 
where = cos^ and = sin^. The diagonalizing of the 4x4 neutral Higgs mass-squared 



matrix has been analytically carried out under some assumptions in Ref. |T5[ and the results 
reduce to Eq. (||) and (|T0|) in the case of large tan /3. 

The another 2x2 matrix can be similarly deal with. Because the couplings of the third 
physical neutral Higgs boson and neutral Goldstone to down-type quarks and leptons are 
not enhanced for large tan/? case in which we are interested, we do not show the explicit 
results. 

Now, we turn to the discussion of the Higgs-fermion-fermion couplings. After completing 
the transformation from the weak states to the mass states, the couplings of neutral Higgs 
to fermions which are relevant to our analysis are 

igruf 



HUf ■ 
HUf ■■ 



cos j3 
igrrif 



(11) 



2mu; cos P 

where / represents down-type quarks and leptons. And the couplings of the charged Higgs 
bosons to fermions are the same as those in the CP-conservative 2HDM (model H, for 



examples see Ref. [|T6|). This is in contrary with the model III in which the couplings of the 
charged Higgs to fermions are quite different from model II. It is easy to see from Eq. (|1T]) 
that the contributions come from exchanging NHB is proportional to y^Gps^c^mj/ cos^/3, 
so that the constaints due to EDM translate into the constraints on sin2,^tan^/9 (1/ cos (3 ~ 



tan/3 in the large tan/5 limit). According to the analysis in Ref. |T^, we have the constraint 



sin2^|tan/5 < 50 (12) 

from the neutron EDM. And the constraint from the electron EDM is not stronger than Eq. 
(0). It is obvious from Eq. (0) that there is a constraint on C, only if tan/5 > 50 and the 
stringent constraint on tan/5 comes out and is tan/5 < 50 when = tt/A. 

III. FORMULA FOR B Xst+t' 

Inclusive decay rates of heavy hadrons can be calculated in heavy quark effective theory 
(HQET) [|l^] and it has been shown that the leading terms in l/mq expansion turn out to 
be the decay of a free (heavy) quark and corrections stem from the order l/rriQ |jT^. In 
what follows we shall calculate the leading term. The transition rate for b — > st^t^ can be 
computed in the framework of the QCD corrected effective weak hamiltonian, obtained by 
integrating out the top quark, Higgs bosons and W^, Z bosons 



H^fj = -^v^,v:xY.c^{^im^^) + Y.GQX^AQ^{^A) (13) 

V-^ 4 = 1 i = l 

where Oi{i = 1, ■ ■ ■ , 10) is the same as that given in the ref. 0, Qj's come from exchang- 



ing the neutral Higgs bosons and are defined in Ref. [|T0|. The explicit expressions of the 
operators governing B XgT^r^ are given as follows: 

O7 = (e/167r2)m,(sL«a^^6fi,)F^,, 
Og = (e/167r2)(sL„7''6L„)f7^r, 

Qi = (eVl67r2)(si„&i?„)(rr), 

Q2 = (eVl6vr')(sL«&/?a)(r75r). (14) 



At the renormalization point /z = mw the coefficients Cj's in the effective hamiltonian 
have been given in the ref. and Cq-'s are (neglecting the 0{tgP) term) 

Zsm^'dw i=HuH2 

z,;)t/t (7^1/ i=Hi,H2 « 
2 

Cq3("^w^) = ^(CqiI"^!^) + C'q2("^m/)), 

2 

CQXmw) = 0, z = 5,---,10 (15) 

where 

= c^, Dh2 = iss:, 



-D/fi — — 2 — ' -°^^2— 2 ' 

Xtlnxt XH±lnxH± — xjnxt 
Xt — I Xh± — Xt ' 



J. _ Xtlnxt XH±lnxH± /^gx 

^ (xt - 1)(xh± - 1) {xh± - Xt){xH± - 1) 



with Xi = mf/ml^. In Eq. (p!5|), are given by 

= ^(-^eci + C5C2), 

= ^(c^Ci + S^C2), 



ci = -xh± + c^xh^c^ + «sg) + s^xh2{s(^ - ic^)-, 

C2 = i {-xh± + s^a;/^^ (s^ - ic^) + c^x^a (c^ + is^)) . (17) 



Neglecting the strange quark mass, the effective hamiltonian (|T^) leads to the following 
matrix element for h st^t~ 

+ 2C^m^,SLia^"''^haf-i^'T + CQ.SLhRfr + CQ^SLhRf-i^r], (18) 



where ||,|j2g 



(19) 



with s = q'^/ml, q = (pr+ +Pt-)^- In (f9|) 9{—,s) arises from the one-loop matrix element 
of the four-quark operators and can be found in Refs. The second term in braces in 

(p!9|) estimates the long-distance contribution from the intermediate, ip" ... In our 

numerical calculations, we choose A;(3Ci + C2) = —0.875 |]22[. 

The QCD corrections to coefficients Ci and Cq^ can be incooperated in the standard way 
by using the renormalization group equations. Although the Ci at the scale /i = 0{mb) have 
been given in the next-to-leading order approximation (NLO) and without including mixing 
with Qi, we use the values of Ci only in the leading order approximation (LO) since no Cq. 
have been calculated in NLO. The Ci and Cq^ with LO QCD corrections have been given 



in Ref. 0. 



Cj{mw) 



. 58 
135 



(^10/23 _ 1) 



29 
189 



{r^^'/^' - l)]C,{mw) 



-^.Ql2CQ,{mw)] 



(20) 



C,{m,) = Cs{m^) + - V-''^'') + " V-''^'')]C2{mw), (21) 

as{mw) 33 87 

Cgim,) = Cgimw), (22) 

CqArnb) = V^^^^^^'CQ^imw), ^ = 1, 2, (23) 



where 7q = —4 is the anomalous dimension of slBr, f3o = 11 — 2nj/3, and 77 = 
After a straightforward calculation, we obtain the invariant dilepton mass distribution 

iir 



dr(5 ^ XsT+T- 



ds 



B(B XrW 



a 



4:t\,,JVthV, 



*|2 



9f2 9+2 9 

Dis) = \CI"\\1 + —){! + 2s) + 4\Cr\\l + —){! + -] 
s s s 



o^2 0^2 

+ \C9Wl + 2s) + —(1 - 4s)] + 12Re(C7C|^^*)(l + — ) 
s s 

+1\CqA\s - At') + 1\Cq,\'s + 6Re{C,C*Qjt (24) 

where s=g^/m^, t=m^/mfe, B{B XJu) is the branching ratio, / is the phase-space factor 
and f(x) = l - 8x' + 8x^ - - 24x^ In x. 

The CP asymmetry for the B — >■ Xgl^l^ and B — > Xsl^l~ is defined as 

We also give the forward-backward asymmetry 

where z = cos 6 and 6 is the angle between the momentum of the B-meson and that of 
in the center of mass frame of the dileptons r+r^. Here, 

E{s) = Re{Cl^^C;s + 2C7Q + Cl"C*Q^t + 2CrCy). (27) 

The CP asymmetry in the forward-backward asymmetry for B X^r+r" and B XsT^t~ 
is defined as 

A' (s)-^^^l^^ (28) 

A(.) + A(s)- ^^^^ 

It is easy to see from Eq. (^) that the CP asymmetry is very small because the weak 

phase difference in C^jCl^^ arises from the small mixing of O7 with (see Eq. (pO|)). In 

contrast with it, A^p can reach a large value when tan/? is large, as can be seen from Eq. 

(p7|) and (p!5|). Therefore, we propose to measure A^p in order to search for new CP violation 



sources. 

Let us now discuss the lepton polarization effects. We define three orthogonal unit 
vectors: 

Pi 

PsX Pi 

Cat 



\Ps X pi I 
er = eV X 

8 



where pi and ps are the three momenta of the i~ lepton and the s quark, respectively, in the 
center of mass of the i~ system. The differential decay rate for any given spin direction 
n of the £~ lepton, where ri is a unit vector in the i' lepton rest frame, can be written as 

dT(n) 1 fdT" 



ds 2\ds 



l + {PLeL + PNeN + PTeT)-n , (29) 



'0 



where the subscript "0" corresponds to the unpolarized case, and Pl, Pt, and P/v, which 
correspond to the longitudinal, transverse and normal projections of the lepton spin, respec- 
tively, are functions of s. From Eq. (^), one has 



dT ^ dV 

m = e,- — in = — e,- 



— (n = Cj) + — (n = -Ci) 
ds ds 

The calculations for the Pj's (i = L, T, N) lead to the following results: 



4sV2^ s ' D{s) ' 
_ 3nt Dt{s) 

- ~2s^^^' ^^^^ 

where 

Dl{s) = Re (2(1 + 2s)Cl^fc; + I2C7C* - QICq^C; - 3sCq,C^J , 

Dm{s) = im [2sCq,c; + sCQ.c'Jf* + sCq.c; + Atc^c; + 2tsCl^f *c^) , 

Dt{s) = Re (^-2C7C; + 4Cl"c; + ^ICrl' - Cg^^C* 

Pi (i=L, T, N) have been given in the ref. [9], where there are some errors in P^ and 
they gave only two terms in D^, the numerator of P/v- We remind that P/v is the CP- 
violating projection of the lepton spin onto the normal of the decay plane. Because P/v in 
B Xgl^l^ comes from both the quark and lepton sectors, purely hadronic and leptonic 
CP-violating observables, such as dn or d^, do not necessarily strongly constrain P/v [0]. So 

9 



it is advantageous to use Pn to investigate CP violation effects in some extensions of SM 
p5| . In the model IV 2HDM, as pointed out above, dn and constrain ^\ sin 2^| tan (3 and 
consequently P/v through Cq. {i = 1,2) (see Eq. (p2D). 



IV. NUMERICAL RESULTS 

The following parameters have been used in the numerical calculations: 

rrit = 175Gev, nib = 5.0Gev, rric = l.GGev, rrir = l.77Gev, r] = 1.724, 

niHi = lOOGev, 112^2 = = 200Gev. 

Numerical results are shown in Figs. 1-9. From Figs. 1 and 2, we can see that the 
contributions of NHB to the differential branching ratio dT / ds are significant when tan (3 is 
not smaller than 30 and the masses of NHB are in the reasonable region, and the forward- 
backward asymmetry A{s) is more sensitive to tan/5 than dT/ds, which is similar to the 



case of the normal 2HDM without CP violation [|l0 . 

The direct CP violation A^p {i = 1, 2) and CP-violating polarization P/v of P — XsT~^t~ 
are presented in Figs. 3-7, respectively. As expected, A^^p is about 0.1% and hard to be 
measured. However, A'^p can reach about 10%. A'^p is strongly dependent of the CP 
violation phase C, and comes mainly from exchanging NHBs as expected. From Figs. 6 and 
7, one can see that P/v is also strongly dependent of the CP violation phase ^ and can be as 
large as 5% for some values of ^, which should be within the luminosity reach of coming B 
factories, and comes mainly from NHB contributions in the most of range of ^. 

Figs. 8 and 9 show the longitudinal and transverse polarizations respectively. It is 
obviously that the contributions of NHB can change the polarization greatly, especially 
when tan (3 is large, and the dependence of Pl on CP violation phase ^ is not significant in 
the most of range of ^. The longitudinal polarization of P — XsT+t^ has been calculated in 
SM and several new physics scenarios ^j. Switching off the NHB contributions, our results 
are in agreement with those in Ref. 0. 

10 



In summary, we have calculated the differential braching ratio, back-forward asymmetry, 
lepton polarizations and some CP violated observables for B X^t^t" in the model 
IV 2HDM. As the main features of the model, NHB play an important role in inducing 
CP violations, in particular, for large tan/5. We propose to measure v4^p, the direct CP 
asymmetry in back-forward asymmetry, in stead of A^p, the usual direct CP violation in 
branching ratio, because the former could be observed if tan/5 is large enough (say, > 30) 
and the latter is too small to be observed. It is possible to discriminate the model IV from 
the other 2IIDMs by measuring the CP- violated observables such as A^p, Pn if the nature 
chooses large tan /3. 
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FIG. 2. Backward-forward asymmetry as function of s, where ^ = 7r/4, solid and dashed lines 
represent tan /3 = 10 and 30, dot-dashed line represents the case of switching off Cg. contributions. 



15 



s 10 

0,2 



-2 



0.15 



0.1 



0.05 



-0.05 



-0.1 



-0.15 



-0.2 




FIG. 3. Aqp as function of ^, where s = 0.8, solid and dashed lines represent tan/3 = 10 and 30. 
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FIG. 4. Aqp as function of ^, where s = 0.8, solid and dashed lines represent tan/3 = 10 and 30. 
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FIG. 5. A'f-jp as function of s, where ^ = 7r/4, solid and dashed hnes represent tan/3 = 10 and 
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FIG. 7. Pn as function of ^, where s = 0.8, solid and dashed lines represent tan /? = 10 and 30, 
dot-dashed line represents the case of switching off Cq. contributions. 
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FIG. 8. Pl as function of ^, where s = 0.8, solid and dashed hnes represent tan/3 = 10 and 
dot-dashed line represents the case of switching off Cq. contributions. 
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FIG. 9. Pt as function of ^, where s = 0.8, solid and dashed hnes represent tan/3 = 10 and 
dot-dashed hne represents the case of switching off Cq. contributions. 
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